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The vibrational density of states (VDOS) of nanoclusters and nanocrystalline materials are derived 
from molecular-dynamics simulations using empirical tight-binding potentials. The results show 
that the VDOS inside nanoclusters can be understood as that of the corresponding bulk system 
compressed by the capillary pressure. At the surface of the nanoparticles the VDOS exhibits a 
strong enhancement at low energies and shows structures similar to that found near flat crystalline 
surfaces. For the nanocrystalline materials an increased VDOS is found at high and low phonon 
energies, in agreement with experimental findings. The individual VDOS contributions from the 
grain centers, grain boundaries, and internal surfaces show that, in the nanocrystalline materials, 
the VDOS enhancements are mainly caused by the grain-boundary contributions and that surface 
atoms play only a minor role. Although capillary pressures are also present inside the grains of 
nanocrystalline materials, their effect on the VDOS is different than in the cluster case which is 
probably due to the inter-grain coupling of the modes via the grain-boundaries. 



PACS numbers: 63.22.-|-m,61.46.-|-w,81.07.Bc,02.70.Ns 

I. INTRODUCTION 

The current interest in nanotechnology makes nan- 
oclusters as well as nanocrystalline materials the sub- 
ject of many research activities. Nanoclusters, i.e., par- 
ticles with diameters in the nanometer range contain- 
ing between a few hundreds and hundreds of thousands 
of atoms, are the building blocks of so-called cluster- 
assembled materials. An understanding of the clusters 
is therefore necessary in order to control manufacturing 
and properties of the resulting materials. Nanocrystalline 
materials are polycrystalline materials with grain sizes 
below 100 nm. The interest in these materials is moti- 
vated by the fact that the small grain sizes give nanocrys- 
talline materials unique physical properties which are dif- 
ferent from those of their coarser-grained counterparts. A 
general introduction to both, nanoclusters and nanocrys- 
talline materials can be found in Ref. 1. 

One possible method for the synthesis of nanocrys- 
talline materials is the inert gas condensation technique. 
In this process, nanoparticles produced from the gas 
phase are deposited on a substrate and subsequently m- 
situ compacted and sinteredi^ This technique as well as 
the similarity of the dimensions of nanoclusters and the 
grains in nanocrystalline materials provide a kind of link 
between both systems and several authors have compared 
vibrational properties calculated theoretically for clusters 
with experimental results of nanocrystalline materials.'ii 
In this work, we present results for the vibrational den- 
sity of states (VDOS) in nanoclusters and nanocrystalline 
materials obtained from molecular-dynamics simulations 
using the empirical tight-binding potentials of Cleri and 
Rosatoi^ The fact, that we have derived the VDOS for 
both types of systems with the same model potential. 



enables us to demonstrate that it is impossible to de- 
scribe the VDOS of nanocrystalline materials in terms of 
the VDOS of the clusters since the dominating effect are 
different in the two kinds of systems. While the VDOS 
of nanoparticles is governed by the large number of sur- 
face atoms and the presence of a rather homogeneous 
compressive capillary pressure, the VDOS of nanocrys- 
talline materials is dominated by large volume propor- 
tion of grain-boundary atoms and the coupling of the 
individual grains via the grain boundaries. 

A crucial point for the simulation of nanocrystalline 
materials is the elaboration of an appropriate model 
with a realistic grain-boundary structure. The most 
widely used approaches to generate models of nanocrys- 
talline materials are the controlled solidification of a liq- 
uid used by Phillpot, Wolf and GleiteriSii and the Voronoi 
construction employed by Chen,- Schi0tz, Di ToUa and 
Jacobsen^SiiS as well as Van Swygenhoven and Caro>ii*i^ 
However, both of these methods give fully dense systems 
with negligible amounts of free volume, whereas the ex- 
perimentally observed densities vary between 75 and 97 % 
of the density of their coarser grained counterpartsii^ii^ 
For this reason we have employed another approach in the 
preparation of our nanocrystalline model systems. The 
models used for this work were obtained from simula- 
tions of the consolidation of nanoparticles into nanocrys- 
talline materials. This approach mimics the experimental 
preparation of nanocrystalline materials by inert gas con- 
densation and leads to porous systems with realistic den- 
sities. The reason why this method is not widely used in 
the literature is that the simulation of the consolidation 
is computationally very demanding. Early applications 
of this techniquei^ii^ were therefore strongly limited in 
the size of the simulated systems. However, the rapidly 
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growing capabilities of modern parallel computers made 
it possible for us to simulate the consolidation of model 
systems with more than one million atoms. 

Experimentally, it has been found that nanocrystalline 
metals exhibit an anomalous increase of the VDOS at low 
energies jiii2ii2i22i2i Similar enhancements have been ob- 
served in computer simulations. ^^■^'^ In Ref. I23L a nickel 
model based on the Voronoi construction has been com- 
pared with an early version of our model for nanocrys- 
talline copper. From this comparison of a fully dense 
model with a porous model, the dominating effect of the 
grain-boundary atoms on the VDOS of nanocrystalline 
metals was deduced. We obtain the same results from a 
comparison of the VDOS of one model system derived at 
different stages of the densification. 

For a discussion of the vibrational properties of 
nanoparticles it is useful to distinguish between the sur- 
face and the core of the particles. In their work on the 
vibrational properties of clusters, Kara and Rahman re- 
ported that the VDOS of the inner atoms of a cluster is 
very similar to the VDOS of the corresponding bulk ma- 
terial but shifted to higher energiesi4 Later, it has been 
shown that this shift is caused by the capillary pressure 
acting inside nanoparticles 1^ Naturally, the distinction 
of surface and core contributions requires a suitable def- 
inition of these regions. In Ref. 'J, the core (or inner 
part) of the cluster was defined to consist of all atoms 
with the bulk coordination number Z — 12. The prob- 
lem with this definition is that a lot of the atoms with 
Z = 12 near the surface show a very different behav- 
ior than those deep inside the cluster. For this reason 
it seems more reasonable to use a definition that gives 
smaller and more homogeneous particle cores. This was 
done in Ref. |3 by considering the second invariant J2 
of the local stress tensor, which turns out to be fairly 
constant inside the particles. The disadvantage of this 
method, however, is that it contains some arbitrariness 
since the onset of the deviation of J2 has to be determined 
by manual inspection. Here, we propose another crite- 
rion based on the common neighbor analysis (CNA)^^ 
which yiel ds very similar results than the method used 
in Ref. |2J without arbitrariness. 

For a macroscopic spherical particle of radius R it is 
known that inside the sphere acts a capillary pressure p 
of magnitude 
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where r is the (average) surface stress of the material. 
Previously, it has been shown for two different types of 
embedded-atom method potentials that this equation re- 
mains valid for particles in the nanometer regimeM^ In 
the present work we find that, for the empirical tight- 
binding potentials of Cleri and RosatO)^ Eq. (Q) holds 
for particles with diameters above a critical size of ap- 
proximately 2.5 nm. Moreover, we report the values of 
the surface stress r for Ag, Au, Cu, Ni, and Ft obtained 
from our simulations with the help of Eq. (^) . 



II. DETAILS OF THE CALCULATIONS 
A. General 

In this article we report on results of molecular- 
dynamics simulations employing the empirical tight- 
binding second moment potentials of Cleri and Rosato^ 
for the calculation of the interatomic forces. These 
potentials, which are fitted to the experimental values 
of cohesive energy, lattice parameters and elastic con- 
stants, have been successfully applied in many works con- 
cerning clustera^''i^^i^^i'^°i^-^i'^^i^'^i^'^ and nanocrystalline 
materials. ii22i2S*^Si2LS The equations of motion in our 
simulations were integrated with the help of the velocity 
form of the Verlet algorithmS^ using a time-step ft, = 2 fs. 
All simulations were carried out at a room temperature 
(T — 300 K). For systems containing more than a few 
thousands of atoms it becomes unfeasible to calculate 
the frequencies of the vibrational modes by diagonaliza- 
tion of the dynamical matrix. For this reason we derived 
phonon density of states p{h(jj) from the Fourier trans- 
form of the velocity autocorrelation function (v(i)v(O)) 
obtained from the simulations using the methods de- 
scribed in the book of Allen and Tildesleyi^ The knowl- 
edge of the VDOS makes it possible to calculate a num- 
ber of thermodynamic properties within the harmonic 
approximations. Here, we use the VDOS to derive the 
specific heat per atom c of the simulated systems with 
the help of 
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B. Nanoparticles 

For the investigation of nanoparticles we have simu- 
lated Ag, Au, Cu, Ni, and Pt nanoparticles with diam- 
eters in the range of 2 to 10 nm. As in our previous 
work,^^ the simulations were started with spherical con- 
figurations of 531, 791, 1205, 1865, 3043, 8247, 17957, 
and 36417 atoms on an fee lattice. During the initializa- 
tion, care was taken to achieve vanishing total momenta 
and angular momenta of the systems. These configura- 
tions were then equilibrated at T = 300 K over a period of 
50 000 simulation steps before the main simulation runs 
of another 50 000 steps. 

After the simulations we used the CNA to group the 
atoms in the final configurations into four classes: per- 
fect fee (FFCC) atoms with an fee environment up to the 
fourth neighbor shell, good fee (GFCC) atoms with an 
fee environment in the first neighbor shell, surface atoms 
(SURF) with coordination number Z < 10, and other 
atoms (OTHR). It turns out that for the 531 and 791 
atom clusters the FFCC group is identical to the core 
found with the previously-used methodi^ Motivated by 
this finding we identify in what follows the particle cores 
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with the PFCC atoms. The advantages of this defini- 
tion are that, on the one hand, it leads to much more 
homogeneous cores than the simple coordination number 
criterion of Ref. '3 while, on the other hand, it avoids the 
arbitrariness of the method employed in Rcf. 24. 

The core definition explained in the preceding para- 
graph enables us to calculate the capillary pressure p act- 
ing inside the particles as a sum over the local pressures 
at the core atoms: 

^'=^Ee(-^«)^+E^>s)- (3) 



TABLE I: Composition of the starting configurations for the 
nanoparticle sintering simulations. Each line specifies the 
number n of clusters containing atoms and their diame- 
ters d in nm for Cu, Ni, and Ag. The total number of atoms 
is 1 007 120 in all three cases. 



n N dcu rfNi dAg 

20 3043 4.0 3.9 4.5 

60 8247 5.6 5.5 6.3 

15 17957 7.3 7.2 8.3 

5 36417 9.3 9.1 10.5 



Here, mi is the mass of particle i while vf, -Fj", and rfj 
are the Cartesian components (a € x,y,z) of the veloc- 
ity of particle i, the force acting between particles i and 
j and their distance, respectively. Furthermore, angular 
brackets denote thermal averaging and the primed sums 
run over all core particles, i.e., all PFCC particles. Fi- 
nally, Nc is the number of particles in the core and Qc 
the average atomic volume of the core particles which we 
have derived from the mean nearest neighbor distance 
fnn inside the core as flc = ^^/2f^^. 

In order to compare the capillary pressures obtained 
from Eq. Q with the prediction of Eq. ^ the radius 
R of the particle has to be known. Since, to our knowl- 
edge, there is no exact and unique definition of the radius 
of a cluster, we adopted here the previously used simple 
scheme for this purposeiSi Thus, the cluster radius R is 
defined to be equal to the largest distance of a particle 
with coordination number Z > 8 from the particles cen- 
ter of mass. 

Finally, we have derived the partial contributions of 
the four groups of atoms (PFCC, GFCC, OTHR, and 
SURF) to the total VDOS of the Cuygi particle. This 
was done by averaging the individual VDOS contribu- 
tions obtained from the velocity autocorrelation function 
of 1000 simulation runs over a period of 25 000 steps, 
each. Such a long simulation time was necessary to com- 
pensate for the poor statistics caused by the small num- 
ber of particles in the system. 



C. Nanocrystalline Materials 

As stated in the introduction the nanocrystalline 
model configurations used in this work were obtained 
from simulations of pressure-assisted nanoparticle sinter- 
ing. With this method we have constructed models for 
nanocrystalline Cu, Ni, and Ag. 

The starting point of the sintering simulations were 
configurations composed of the equilibrated configura- 
tions of clusters described in the preceding section. Ran- 
domly oriented copies of these clusters were placed in a 
cubic box with twice the volume of a bulk system with 
the same number of atoms. In order to avoid overlapping 
of the clusters in the initial configurations, the positions 



of the clusters were determined with the help of a Monte- 
Carlo procedure. In this procedure, the clusters were re- 
garded as point particles interacting by a purely repulsive 
pair potential. *° For each pair of particles the range of 
the potential was chosen to be larger than the sum of the 
particle radii plus the interaction range of the Cleri and 
Rosato potential used in the molecular-dynamics simula- 
tions. Due to this construction the interaction energy in 
the Monte-Carlo simulation is exactly zero if all clusters 
are well separated and positive otherwise. Starting from 
arbitrary positions the clusters were then moved in trial 
steps which were only accepted if they resulted in a lower 
total interaction energy. This procedure was continued 
until the total energy remained zero for 100 continuous 
trial moves of all clusters. Details of the numbers and 
sizes of the clusters used for the generation of the initial 
configuration are given in Table 

After their generation the starting configurations were 
simulated at constant temperature and pressure using the 
ParincUo- Rahman scheme^ restricted to orthogonal sim- 
ulation boxes and the Nose-Hoover thermostat methodi^ 
The application of the latter method was necessary to 
remove the excess energy released during the densifica- 
tion from the system. The simulations were started with 
a simulation run at zero pressure over 50 000 simulation 
steps (100 ps). Another 50 000 simulation steps were then 
performed at a hydrostatic pressure of 100 MPa. This 
predensification phase was followed by three consecutive 
simulation runs at pressures of 0.7, 1.4, and 2.1 GPa with 
a length of 100000 steps (200 ps), each. In the last step, 
the intermediate configurations resulting from the sim- 
ulations at 0.7, 1.4, and 2.1 GPa were slowly relaxed to 
zero pressure, again. This was done in a simulation where 
the pressure was reduced by 175 MPa after each 2500 sim- 
ulation steps. Eventually, the simulations of the relaxed 
systems were continued for 10 000 steps at zero pres- 
sure in order to diminish possible perturbations caused 
by the stepwise reduction of the pressure. The nominal 
grain sizes of the final configurations, calculated under 
the optimistic assumption that each of the initial parti- 
cles forms one grain in the resulting material, are 7nm in 
the case of Ag and 6 nm for Cu and Ni. The actual val- 
ues, however, will be somewhat higher since during the 
compaction surely some of the particles joined to form 



4 



single grains. 

The models resulting from the relaxation of the the sys- 
tems after the sintering under different pressures enabled 
us to analyze the VDOS of the nanocrystalline metals at 
different densities. We did this analysis for the models 
of nanocrystalline copper obtained from maximum pres- 
sures of 0.7 and 2.1 GPa. For these configurations we de- 
rived the VDOS from the velocity autocorrelation func- 
tion averaged over 10 consecutive simulation runs with 
a length of 10 000 steps, each. In addition to the total 
VDOS, we calculated the partial contributions of the fol- 
lowing groups of atoms, identified by a CNA analysis of 
the averaged atomic positions: atoms with a perfect fee 
or hep environment up to at least the fourth neighbor 
shell (PFCC, PHCP), good fee or hep atoms with an fee 
or hep like environment in the first and more neighbor 
shells (GFCC, CHOP), surface atoms with coordination 
number Z < 10 (SURF) and grain boundary atoms (GB) 
not belonging to any of the other groups. 

Nanocrystalline materials inherently contain a certain 
amount of excess volume in the form of open or closed 
pores. In order to analyze the free volume in our model 
configurations we followed a procedure similar to that de- 
scribed by Campbell et al.4i we subdivided the simula- 
tion boxes into a regular grid of voxels (the 3-d equivalent 
of a pixel) and performed an analysis of the clusters of 
empty voxels. However, in contrast to Campbell et al., we 
did not use simple cubic voxels for this analysis but the 
Wigner-Seitz cell of the fee lattice since this choice gives 
considerably less anisotropic voxels. The edge length of 
the voxels were chosen to be 5.75, 5.60, and 6.50 a. u. for 
Cu, Ni, and Ag, respectively. Incomplete voxels, result- 
ing from the incommensurability of the simulation boxes 
with respect to the voxel lattice, were excluded from this 
analysis. 



III. RESULTS 
A. Nanoparticles 

In Fig. ^ we present the capillary pressures inside the 
nanoparticle configurations calculated from our simula- 
tions with the help of Eq. and fitted to the simple 
model of Eq. . From this figure it is clear that Eq. 
accounts for the capillary pressure in all configurations 
except the smallest ones. Thus, within the Cleri and 
Rosato potentials the scaling regime where the capillary 
pressure scales monotonously with the system size ex- 
tends down to particle sizes of 2 to 3 nm. For the small- 
est clusters below this size, Fig.^shows strong deviations 
from the linear behavior of Eq. ^ whose signs and mag- 
nitudes vary strongly between the elements. 

The fact that the capillary pressure inside the sim- 
ulated nanoparticles follows the linear relationship of 
Eq. Q makes it possible to derive the average surface 
stress T for the different elements. The results are pre- 
sented in Table Altogether with some experimental data. 
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FIG. 1: Capillary pressure inside fee metal nanoclusters as 
a function of the inverse cluster radius BT^ . The lines are 
the results of linear fits to the data of all but the smallest 
particles. 



Two points have to be noted here. First, in the case of 
Cu, there is a remarkable difference between our result 
and the experimental value. However, compared to the 
other metals, it is the experimental value which seems to 
be exceptionally low. In fact, in Ref. 23 two different 
surface stresses of 0.0 and 5.3 N/m have been reported. 
These values follow from lattice constant measurements 
using two different electron diffraction spots. The au- 
thors of Ref. 113 reject the higher value since the corre- 
sponding observations of the (111) peak were less reliable 
than those of the (220) peak which gave the lower value. 
The second point to be noted is the strong scatter in the 
experimental data of Au and Pt. One possible reason for 
this are the relatively large error bounds of up to 1 N/m. 
In addition to this, Solliard and Flueli report that their 
values for the surface stress of Au have been increased 
by the occurrence of structural changes of the clusters4& 
However, consistent calculation of the surface stress from 
Eq. (Q) is only possible if all data points correspond to 
the same structure. Taking all these experimental uncer- 
tainties into account, the surface stresses resulting from 
our simulations are in good agreement with the experi- 
mental data. In order to shed more light on the case of 
Cu, independent measurements or ab-initio calculations 



TABLE II: Calculated and experimental values of the surface 
stress r in various fee metals. 
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FIG. 2: Normalized VDOS of the CU791 cluster, (a) Total 
VDOS compared to the surface (SURF) and core (PFCC) 
contributions; (b) Core part of the VDOS compared to the 
VDOS of crystalline fee Cu under pressures of 0. and 1.5 GPa; 
(c) Comparison of the three non-core VDOS contributions. 



would be highly desirable. 

The VDOS of the Cuygi particle and its partial contri- 
butions are shown in Fig.|3 In the upper panel, Fig.l^Ja), 
the contributions from surface and core atoms are com- 
pared to the total VDOS, demonstrating the huge differ- 
ences between both contributions. While the core contri- 
bution has two clearly separated peaks corresponding to 
the longitudinal and transversal modes of the fee VDOS, 
the surface contribution has only one peak at energies 
below the transversal peak of the core and a shoulder at 
the position of the longitudinal peak. The form of the 
surface VDOS resembles closely the VDOS found near 
flat surfaces.'*^ From these results, it is clear that the 
surface atoms lead to a significant increase of the VDOS 



of nanoparticles at low energies whereas they reduce the 
VDOS at the upper end of the spectrum. 

In Fig.|2Ib) the particle's core VDOS is given together 
with the VDOS of crystalline fee Cu under pressures of 0. 
and 1.5 GPa, where 1.5 GPa corresponds to the capillary 
pressure acting on the cluster's core. This comparison 
reveals, in accordance with previous findings, that the 
VDOS inside the core of a nanoparticle is well accounted 
for by the bulk system under pressure. In particular, the 
capillary pressure gives rise to a small shift of the longitu- 
dinal peak to higher energies, thereby slightly increasing 
the high-energy VDOS. 

Finally, Fig. ^c) shows the three non-core contribu- 
tions to the VDOS of the CU791 particle. From this fig- 
ure, it can be seen that there is a gradual transition from 
the GFCC to the SURF contribution which reflects the 
geometry of the three groups of atoms: The perfectly fee 
ordered core is directly surrounded by a shell of GFCC 
atoms. The GFCC atoms in turn are surrounded by 
subsurface atoms with coordination numbers Z — 10, 11 
followed by the outermost layer of surface atoms with 
Z < 10. 



B. Nanocrystalline Materials 

During the simulations of nanoparticle sintering, the 
systems evolved from a gas of individual particles via a 
loose agglomerate of particles to a dense solid material. 
An important parameter of this compaction process is 
the relative density defined as the actual density of the 
system divided by the density of the corresponding crys- 
talline bulk material. From the relative density it can be 
seen how much excess volume remains in the compacted 
system. In Table IIIII we report the relative densities of 
our simulated nanocrystalline materials after the sinter- 
ing pressure had been removed. Not surprisingly, the 
relative density increases for all three elements as the 
sintering pressure (and time) increases. However, there 
are some subtle differences between the metals which are 
worth noting. After compaction without application of 
an external pressure Cu and Ni both show the same rel- 
ative density of 77% whereas the Ag system has a con- 
siderable lower density of 70%. In contrast to this, the 
picture has changed at the end of our simulations, where 
the Ag system has the highest density, followed by Cu 
and Ni. 

The different behavior of the three metals during the 
compaction described in the previous paragraph can be 
understood from the differences of the stacking-fault and 
surface energies of the three metals. At the beginning of 
the compaction, when the system is not a solid system 
but a loose powder, the compaction rate is mainly deter- 
mined by the surface energy of the metals which is the 
driving force behind the whole sintering process. Since 
the surface energy is highest for nickel (see e. g. Ref. 
and lowest for Ag, the densification starts fastest for Ni 
and slowest for silver. However, as the compaction pro- 
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ceeds, further densification requires the plastic deforma- 
tion of the particles which is controlled by the stacking- 
fault energy. The densities at the end of the simulations 
therefore reflect the fact that (experimentally as well as in 
the model) the stacking- fault energy in Ag is lower than 
in Cu which in turn is much lower than that in Ni42i The 
influence of the stacking-fault energy on the compaction 
rate can also be seen from the number of atoms with an 
hep environment, which indicate the presence of stacking 
faults due to plastic deformation. Among the configura- 
tions belonging to the first line of Tablc lTTll the Ag system 
has the lowest number of hep atoms. This means that 
the system has undergone less plastic deformation than 
the other configurations and explains the rather low rela- 
tive density of the silver system at this point. At the end 
of the simulations, however, the situation has reversed. 
Due to its high stacking-fault energy, it is now the Ni 
configuration which has the lowest number of hep atoms 
as well as the lowest relative density. 

While the relative density gives information about the 
amount of excess volume remaining in the systems, it 
tells nothing about its structure or distribution. Such 
information can be gained from the free volume analy- 
sis described at the end of section III CI Figure 13 shows 
the free volume found by this method in our models of 
nanocrystalline copper after sintering under pressures of 
0.7 and 2.1 GPa with relative densities of 89 and 96%, 
respectively. It can be seen that there is a significant dif- 
ference between the pore structures of both systems. In 
the less dense system, most of the free volume is stored in 
one big percolating cluster, whereas in the denser system 
the free volume is distributed over several isolated clus- 
ters. This means that during the sintering the structure 
of the free volume changed from an extended network of 
open pores to a closed pore state. Analysis of all config- 
urations shows the occurrence of extended open pores in 
all configurations with relative densities up to 90 % and 
closed pores at higher densities. 

In Fig. ^ the total VDOS of the two Cu systems dis- 
cussed in the previous paragraph are compared with that 
of crystalline fee copper. The VDOS of both nanocrys- 
talline systems show strong enhancements at high and 
low energies, as seen before theoreticallj*22i as well as 
experimentallyiSS The important point here is that the 
curves of both nanocrystalline systems are nearly identi- 



TABLE III: Relative densities of the nanocrystalline configu- 
rations after removal of the maximum sintering pressure Pmax 
(GPa). 



Pmax Cu Ni Ag 

0.0 0.77 0.77 0.70 

0.7 0.89 0.85 0.90 

1.4 0.94 0.89 0.96 

2.1 0.96 0.92 0.98 




FIG. 3; Free volumes in nanocrystalline Cu configurations 
after sintering with relative densities of 89 % (top) and 96 % 
(bottom). Clusters containing less than three voxels have 
been suppressed. 



cal, although Fig. O shows clearly a big difference in the 
amount of internal surface area in both systems. This 
establishes that the enhancements are not caused by any 
kind of surface effects. 

In addition to the total density of states we also cal- 
culated the partial contributions to the VDOS for our 
densest copper system. The results of these calculations 
are shown in Fig.|Sl In panel (a), we compare the VDOS 
from atoms with perfect (FFCC) or good (GFCC) fee 
environment, i.e., atoms inside the grains, with the crys- 
talline VDOS. It can be seen that the VDOS inside the 
grains of the nc material is very similar to that of the 
crystalline bulk system. The difference consists of a slight 
broadening and a loss of the sharp structure of the peaks 
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FIG. 4: Normalized total VDOS of nanocrystalline (nc) cop- 
per at 89 and 96 % relative density together with that of crys- 
talline (fee) copper. 



due the a lack of long range order and the presence of 
a strong stress field. Interestingly, the peak of the lon- 
gitudinal modes is not at all shifted with respect to the 
crystalline VDOS. This is striking, since we find that a 
net pressure of 460 MPa acts on the PFCC atoms forming 
the core of the grains. Although lower than the capillary 
pressure acting inside the CU791 cluster one would ex- 
pect that this pressure gives rise to a visible shift of the 
longitudinal modes. A possible reason for the observed 
different behavior is that the modes inside the grains are 
extended, i.e., they are coupled among the grains via the 
grain boundaries. In that case the looser couphng be- 
tween the grains might counteract the stronger coupling 
inside the grains due to the capillary pressure. 

Although the broadening of the intra-grain VDOS 
leads to a small increase of the VDOS at high and low 
energies, the effect is not strong enough to account for 
the strong increase of the total VDOS visible in Fig. 0] 
The biggest contribution to this increase stems from the 
high amount of grain-boundary atoms which show a very 
broad and nearly structureless VDOS with strong en- 
hancements at high and low energies (see Fig.jSJb)). In 
contrast to this the VDOS of atoms at internal surfaces, 
also given in this figure, shows an even stronger increase 
of the VDOS at low energies but no increase in the high 
energy regime. This gives further support to the obser- 
vation that the contribution of the surface atoms to the 
observed enhancements of the total VDOS is insignifi- 
cant. 



C. Nanoparticles vs. Nanocrystalline Materials 

Having discussed our results concerning nanoparticles 
and nanocrystalline materials, we now compare directly 
the vibrational properties of the two kinds of systems. 
For this purpose we present in Fig.|H|the total VDOS of 
the Cuygi cluster and nanocrystalline copper sintered at 
2.1 GPa together with the VDOS of crystalline fee cop- 
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FIG. 5: Normalized partial VDOS of nanocrystalline (nc) cop- 
per at 96 % relative density. See text for details. 



per. From this figure it can be seen that both nanoscale 
systems show a broadening of the VDOS which gives rise 
to an increase at high and low phonon energies. However, 
while for the nanocrystalline material the increase is sim- 
ilar on both sides of the spectrum, it is much stronger 
at low than at high energies in the case of the cluster. 
The reason for this difference lies in the nature of the 
atoms contributing most to the increase of the VDOS. As 
discussed above, the low and high energy enhancements 
of the VDOS in the nanocrystalline material are mainly 
due to the high number of grain boundary atoms which 
show a rather symmetric broadening of the VDOS (see 
Fig. Etb)). In contrast to this it is clear from Fig. [21 a) 
that the VDOS enhancements in the cluster are driven 
by two different types of atoms: The strong low energy 
increase reflects the VDOS contribution of the surface 
atoms whereas the high energy increase is caused by the 
shift of the core VDOS due to the capillary pressure. Fi- 
nally, it can be seen from Fig.|Hlthat the shift observed for 
the longitudinal modes in the cluster core only leads also 
to a visible difference between the VDOS of the cluster 
and the nanocrystalline material. 

In order to study the effect of the VDOS changes in 
the nanoscale materials on other physical properties, we 
have calculated the specific heat of our model systems 
with the help of Eq. |j2Jl. The results presented in Fig.[7| 
reveal an increase of the specific heat in the nanoscale 
systems at low temperatures. All three nanoscale sys- 
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Phonon energy (meV) 

FIG. 6: Normalized VDOS of a CU791 cluster, nanocrystalline 
(nc) copper with a relative density of 96 % and crystalline (fee) 
copper. 

terns show the maximum of the excess specific heat with 
respect to the crystalline state at a temperature around 

50 K. The height of this maximum, on the other hand, 
is much stronger in the case of the nanoparticle due to 
the stronger increase of the VDOS at low energies in this 
system. Our results here are in good qualitative and 
quantitative agreement with other theoretical results for 

51 nanoparticles'* and a nanocrystalline Lennard-Jones 
system. Experimentally, excess specific heats have been 
observed for several nanocrystalline metalsi^2*^ii^24^i^ 
At least in some of these works )^2i^ the measured in- 
crease of the specific heat is considerably higher than the 
theoretical results. However, it has been pointed out by 
Tschope and Birringer^ that large part of the measured 
excess specific heat might be due to contaminations of 
the samples with lighter elements. Loffler,^^ on the other 
hand, reports increases of the specific heat by approxi- 
mately 7 % in nanocrystalline Pd, which compares well 
to our results although he finds the maximum of the in- 
crease at a temperature slightly above 150 K. 

IV. SUMMARY AND DISCUSSION 

In this work we have studied nanoparticles and 
nanocrystalline materials made of fee metals with the 
help of molecular-dynamics simulations based on the em- 
pirical potentials of Cleri and Rosatoi^ 

In order to distinguish between core and surface parts 
of the nanoparticles, we made use of the common neigh- 
bor analysis^^ and defined the core of the particles as the 
set of atoms with a perfect fee environment. With the 
help of this criterion we studied the capillary pressure 
acting inside spherical nanoparticles of Ag, Au, Cu, Ni, 
and Pt with diameters between 2 and 10 nm. We find 
that the capillary pressure follows the same behavior as 
macroscopic spherical particles down to sizes of approx- 
imately 2.5 nm. These findings allowed us to derive the 
average surface stress values for Ag, Au, Cu, Ni, and Pt 




■ 50 100 150 200 250 300 
Temperature T{K) 

FIG. 7: (a) Specific heats calculated for crystalline (fee) cop- 
per, nanocrystalline copper with relative densities of 89 and 
96% and a CU791 cluster; (b) Excess specific heat of the 
nanoscale systems in panel (a) relative to the crystalline state. 

predicted by the potentials of Cleri and Rosato. Except 
for Cu our results agree with the available experimental 
data. 

Our investigations of the VDOS of a Cuygi nanopar- 
ticle show, in agreement with our previous findings;^ 
that the VDOS inside the core of the cluster is identical 
to the VDOS of fee copper under the capillary pressure. 
One particular effect of the capillary pressure is a shift of 
the longitudinal modes peak to higher energies. Outside 
the core the VDOS changes gradually to the typical sur- 
face behavior of the outermost atoms. Compared to the 
VDOS of the bulk material, the cluster VDOS appears 
to be slightly broadened with an increased VDOS at high 
and low energies due to the core and surface atoms, re- 
spectively. 

The model configurations for nanocrystalline Ag, Cu, 
and Ni used in this work have been obtained from exten- 
sive simulations of the sintering of nanoparticles. These 
calculations resulted in nanocrystalline configurations of 
the three metals with different relative densities in the 
range 70 to 98%. The development of the relative den- 
sity of the three metals during the compaction process 
suggests that the compaction rate of the material is ini- 
tially governed by the surface energy of the material. As 
plastic deformation becomes more and more important 
during the compaction, the stacking-fault energy starts 
to dominate the process. It should be noted, however. 
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that the role of other mechanisms Hke surface or grain- 
boundary diffusion might be underestimated by our sim- 
ulations since such processes are beyond the time-scale 
of molecular-dynamics simulations. Further analysis of 
the excess volume remaining in the nanocrystalline con- 
figurations revealed a transition from an extended open 
pore network to closed pores at relative densities above 
90% for all three metals. 

In Ref. a model of nanocrystalline Ni based on the 
Voronoi construction has been compared with an early 
version of our model of nanocrystalline copper. This 
comparison showed that the experimentally observed en- 
hancement of the VDOS of nanocrystalline materials at 
high and low energies is due to grain-boundary atoms, 
whereas the influence of atoms at internal surfaces is neg- 
ligible. In the present work we obtain the same result, 
using one and the same model of nanocrystalline cop- 
per for all the calculations. Although we find a small 
pressure to be present inside the grains of our model sys- 
tems, we do not observe a shift of the longitudinal modes 
in the nanocrystalline material. A possible explanation 
for this difference between the dense nanocrystalline ma- 
terial and the nanoparticle is the loose coupling of the 
modes inside the grains via the grain boundaries. 

Our calculations of the specific heat of nanocrystalline 
Cu and a Cu nanoparticle show that the low energy 
VDOS enhancements in the nanoscale systems give rise 
to a pronounced peak in the excess specific heat. This 
result is in qualitative and quantitative agreement with 



other theoretical calculations^-^- and at least in qualita- 
tive agreement with experimental findings i^SiSiiSiSSiSi 

This work provides a deeper understanding of the vi- 
brational properties of nanoclusters and especially the 
differences between these two nanoscale systems. For 
future work it would be important to investigate how 
the strong stress field present inside the nanocrystalline 
model systems affects the vibrational properties as well 
as some general study of the properties of crystalline sys- 
tems coupled by a disordered medium. Finally, it would 
be interesting to use some multi-scale method to relax 
our model systems beyond the time scale of molecular- 
dynamics simulations. 
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